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The efTect of an electric field on conduction in a disordered system is an old but largely un- 
solved problem. Experiments cover an wide variety of systems - amorphous/doped semiconductors, 
conducting polymers, organic crystals, manganites, composites, metallic alloys, double perovskites 
- ranging from strongly localized systems to weakly localized ones, from strongly correlated ones 
to weakly correlated ones. Theories have singularly failed to predict any universal trend resulting 
in separate theories for separate systems. Here we discuss a recent one-parameter scaling that has 
been found to give a systematic account of the field-dependent conductance in two strongly localized 
systems of conducting polymers and manganites. The nonlinearity exponent, x associated with the 
scaling was unexpectedly found to possess multiple values. We find that the scaling applies to all 
other systems as mentioned. For 3D strongly localized systems the exponent lies between -1 and 
1, and surprisingly, is quantized {x ^ 0.08 n). For 2D weakly localized systems, the nonlinearity 
exponent a; is ^ 9 and is roughly inversely proportional to the sheet resistance. Existing theories 
of weak localization prove to be adequate in explaining the proposed scaling which however poses a 
challenge for theories in case of strong localization. 

PACS numbers: 66.30.Pa, 62.20.Dc, 81.07.-b 



I. INTRODUCTION 



Effects of disorder on transport phenomena have been 
studied for more than four decadesi, yet understanding 
is still far from complete. One of the effects of disorder is 
to drive samples into nonOhmic regime even under small 
bias. This effect is particularly pronounced in small-sized 
samples such as nanodevices. It is thus necessary to un- 
derstand the phenomena of nonlinear transport to ensure 
that the devices operate as desired vifith optimum prop- 
erties. In seventies the focus was on the materials such as 
amorphous/doped semiconductors which are recognized 
today as strongly localized systems characterized by the 
exponential dependence of conductivity on temperature. 
The conduction occurs through variable range hopping 
(VRH) between randomly occurring localized states. The 
theoriesi"— proved inadequate to account for the exper- 
imental data in the full range of applied field (see be- 
low and Ref. [3] for a recent critique). Subsequently, 
new disordered regimes, namely weak localization in 2D 
and metal-insulator transition (MIT) in 3D were identi- 
fied by Abrahams et ali^. The experiments^ were mostly 
limited to verification of the prediction^ of logarithmic 
variation in weak localization of conductivity with either 
temperature at small bias or large bias at low temper- 
ature. Here we attempt to find a universal description 
of the response to an applied electric field by adopting a 
model-independent scaling approach. The strong motiva- 
tion comes from the remarkable similarity in the response 
of different disordered regimes in diverse materials such 
as mentioned earlier. In this approach the quantity which 
is of major interest is the field scale (or the onset field), Fq 
at which a sample deviates from the Ohmic behavior. It 
is shown that while the scaling formally is same across the 
regimes the associated exponent x varies and is nonuni- 



versal. In 2D weak localized regimes, the exponent x 
is roughly inversely proportional to the sheet resistance. 
It becomes quantized in strong localization. The cur- 
rent theories of weak localization and region around the 
MIT support the scaling. This makes the scaling quite 
plausible in strong localization which however remains a 
theoretical challenge. All discussion here are limited to 
two and three dimensions only. 



II. SCALING OF FIELD-DEPENDENT 
CONDUCTANCE 

The field-dependent conductance can be written in 
general as S](M,F) = T,o{M) <^{F,M) where F is the 
field and M stands for the control parameter (most com- 
monly, temperature) which is varied to change the lin- 
ear conductance So(M) = E(Af, 0). Following the spirit 
of the scaling theory^ we assume that <I> is a function 
of F and T,o{M) only and not of temperature or size 
explicitly. This is also motivated by the goal to have 
quenched disorder (corresponding to samples fabricated 
under different conditions) as a valid parameter for M. 
Thus, S(Af,F) = T.o{M) $(-F, I]o(M)). It is natural to 
expect that $ will be a function of F/Fq. Hence the scal- 
ing that the disordered systems are expected to follow is 
given by^i^ 
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In analogy to power-laws in critical phenomena we have 

F,{M) ^ j:,{Mf" , (2) 

where the 'critical' conductance is zero for all cases, xm 
is the nonlinearity exponent with a subscript for the con- 
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trol parameter. The above relation ahows seamless con- 
sideration of disorder as the control parameter (see Fig. 
6b of Ref. i for an illustration). A lower bound of 
xm may be obtained by considering the current scale 
lo ~ Soi^o So^"*"^. Since Iq should remain finite 
when So is very small this obtains a; ^ — 1. At small 
q[= F/Fo) s: 1, « 1. At large field {q > 1), the 

conductance often becomes independent of T (or, T,o{T)) 
(see Fig. 1) when the control parameter is temperature. 
It is seen from Eq. ^ and ^ that this requirement is 



readily ensured if $((7) 
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at large q ioT xt > 0. 



Thus, at large fields the conductivity varies as a power- 
law, S(r, F) - F""^ such that for xt > 



(3) 



The above relation can be understood by simply writing 
the conductance as E w + bF^^^^ where is a 
constant. This immediately generates Fq ^ '^o^'^ ■ 
Note that the similar scaling (Eqs. 1-3) with F and 
Fo replaced by the frequency and its characteristic 
valued also describes the ac-conductivity at different ac 
amplitudes. 



A. Weak localization 

We first discuss nonlinear conduction in 2D weakly lo- 
calized thin films by revisiting the data of Dolan and 
OsherofS shown in the inset (a) of Fig. 1. The data con- 
tain a lot more information than the simply logarithmic 
variation with the applied field. Indeed, the sample re- 
sistance i?(=l/E) with an initial Ohmic value {Ro) looks 
similar to those in strongly localized systems^. Clearly, 
there exists a voltage scale 14 (T) at each T such that 
R ^ Ro for q{= V/Vo) < 1 and R decreases for q > 1. 
As shown in Fig. 1 (main panel), the collapse of all 
data at various temperatures on a single curve^^ con- 
firms the universality of the scaling in Eq. ([1]). Here 
Vo{T) is chosen such that R{T, Vo) = (1 + e)Ro{T) with 
e = —0.05. The solid line in Fig. 1 indicates a power- 
law at large fields, R ~ V^'^ where the exponent zt is 
1/64.9. The inset (b) displays a log-log plot of Vo vs. Ro 
which yields a large valuoii for the nonlinearity exponent 
XT = 64.9 ±0.4 fulfilling the requirement zt ~ I/xt (Eq. 
[2]). At large fields the resistance variation in weak scat- 
tering regime is given by the hot electron model^i^ 



R(T, V) 
{l + e)Ro{T) 



- 1 = -i?ni;„Ay In 



(4) 



where i?n is the temperature-independent sheet resis- 
tance, E„ = e^/2TT^h = 1.2 X 10-5 r^-i, Ay = ap/{l + 
p/2), and a an universal constant. The exponent p is de- 
fined by Tin ~ where t^ is the inelastic relaxation 
time. The above expression is compatible with Eq. ([T]). 
Considering that the change in R is small, the above can 
be recast as a power-law R/Ro = (1 -t- e){V/Voy^ with 




FIG. 1. Scaling in 2D weakly localized metallic thin films 
of AuPd. Plots of resistance vs. electric bias at different 
temperatures are shown in the inset (a) (data adapted from 
Ref. 13) ■ Scaling of the same data leads to the data collapse 
as shown in the main panel. Ro(T) and VoiT) are the Ohmic 
resistance and the onset bias respectively at temperature T. 
Inset (b) shows log-log plots of Vo vs. Ro (solid symbols) and 
T (open symbols). The solid lines are linear fits to the data 
with the slopes as mentioned. 



the exponent zt = Ru^u^v- With i?n = 4600 Q. and 
Xv = 0.28 (Ref. 0), one determines zt ~ 1/62.9 which is 
close to 1/64.9 obtained from the digitized data (Fig. I). 
To obtain xt from the hot electron mode£ we consider 
the relevant energy balance rclatiooi^: 



(5) 



where is the effective electron temperature and ctq 
is the temperature-independent Ohmic conductivity. At 
the onset of nonlinearity Te T so that 



Fo T^+P/'^. 



(6) 



The above provides the most direct method for deter- 
mining p. A plot of Fo vs. T in the inset (b) yields 
1 +p/2 w 2.4 close to 2.6 reported in Ref. A similar 
expression as in Eq. (|4|) can be written for temperature at 
small bias leading to Ro ^ y-^nS„ap^ Using this in Eq. 
© gives Eq. © with xt = l/Rn^uXy = ^/zt- With 
Ay being of the order of unity, xt ~ 81000/i?n ;» 1 since 
the maximum metallic sheet resistance is predicted to be 
about 30 kO^. Thus, the exponent in the 2D metallic 
state should have the minimum value a;™™ « 3. 

To verify that xt oc l/i?n, exponents (open symbols) 
processed from the data available in literature for four 
different materials (films) of thicknesses 3-13 nm are 
shown in Fig. 2. All data are seen to lie on a line 



3 



10' 



. V 


1 1 

•\ 3-13 nm 


1 ' 










1 






Ax 






140-150 nm~ 


^ - . 





10° 10' 10^ 10^ 10* 10' 
R (CI) 

FIG. 2. EfTect of film tliickness on the nonlinearity exponent 
in 2D metallic states. The exponent is plotted against sheet 
resistance for various materials and film thicknesses. One set 
of data (open symbols) belongs to thin films (D-AuPd [6], A- 
Pt[il, o-Au[3|), and o-Cu[l|]). Another set of data (closed 
symbols) is from thick films (Refs. ITsI ). The solid line is a 
least square fit and the dashed line is drawn parallel to it. The 
vertical line corresponds to the maximum metallic resistance 
of 30 kn. 



with a slope of -0.85 which is close to the expected 
value of -1. The small discrepancy may be due to the 
variation in experimental p. Also shown in the figure 
are the data (closed circles) from two thick filmsi^ 
(Gco.79Auo.21 and CxCui_x of thickness 140-150 nm 
near, but on the metallic side of, the metal-insulator 
transition in 3D^). These thick films appear to be in the 
cross-over regime from 2D to 3D. The reported small 
temperature variation of conductivities (~ 7^0.02-0.07^ 
is akin to the logarithmic one characteristic of a 2D 
system. A film is supposed to be 2D if the inelastic 
mean free path is greater than the thickness. Using 
an alternative of Eq ([5]), eFolrf, ^ kT and Fo ~ 9 V/cm 
at 1.5 K, we estimate ~ 160 nm i.e. of the order 
of film thickness. While these data from thick films 
separately follow the same trend with sheet resistance 
as thin films, the exponents are roughly an order of 
magnitude less than those in thin films for similar 
sheet resistances. This trend in the exponent value 
is in agreement with the experimental results in 3D 
strongly localized systems where xt decreases to a value 
less than 1 (see Tables in Appendix B and below), p 
determined from the onset field data using Eq. ^ is 
0.9 and 0.7 for the two films respectively. These re- 
sults indicate a possible renormalization of the constant 
a with the film thickness. The fitted line yields ~ 9. 



B. Metal-Insulator transition (MIT) 

We now consider a regim o^^'^^ somewhat interme- 
diate between weak localization with the logarithmic 
dependence, and strong localization with the activated 
exponential dependence, of the Ohmic conductivity do. 
Such a regime occurs in the vicinity of the MIT with 



large localization length^. The conduction takes place 
through hopping of electrons across the localization 
centers. The conductivity has a power-law dependence 
on temperatureii, ctq ~ yp/^ in contrast to the ex- 
ponential variation. In do ^ —[Tq/T) . Here To is a 
characteristic temperature and the exponent m can 
vary from 1/4 (variable range hopping or VRH) to 1 
(nearest neighbor hopping)^. Usually, To 3> T. However 
when To ^ T, the average hopping length is of the 
order of the localization length and the exponential 
form is replaced by a power-law. In the case where 
T is the control parameter, we apply the hot electron 
model (Eq [5]) as before and use cto ^ TP/^ to obtain 
ZT = 1/(1 +p/3) and XT = 1 + p/3 = 1/zt (Eq. ^. 
With p ^ 1, XT is predicted to be about 1.3. Rosenbaum 
et alJ^ measured field-dependent conductivity in Si:P 
samples near but on the insulating side of, the MIT 
at different temperatures. The behavior of (JoiT) was 
rather ambiguous. Nevertheless, the data were found 
to obey the scaling (Eq. [1]) to give xt = 0.15 ± 0.01 
(Appendix A). While Uo{T) yielded an unphysically 
large p w 30, Fo{T) yielded p « 1 (Eq. ^ indicating that 
the hot electron model perhaps may not be applicable 
in this regime. Clearly, the experimental exponent (< 1) 
does not match the predicted value (> 1) but the value 
is compatible with the ones found in strong localized 
systems (see below). 



C. Strong localization 

The diverse strongly localized systems as mentioned in 
the beginning are described by the universal scaling as in 
Eq. ID). However a much desired theoretical framework 
which would be independent of the microscopic details as 
in 2D is missing at present. Instead, there exist separate 
theories for separate systems^"— , often falling short of ful- 
filling the general requirements of the model-independent 
scaling. For example, the existing literaturei"— on VRH- 
type transport deals with two field scales. The recent 
stud}*^ has clearly established the fact that there exists 
only one field scale, not two and the scale is not simply 
related to either of the two length scales - the average 
hopping length and the localization length (see also Ref. 
ITsI ). The same is observed in other VRH classes such as 
amorphous/doped semiconductorsi^, composites at low 
temperature^^ whose data are available in literature. 

As a further example of implications of the model- 
independent scaling, let us discuss a non-VRH, corre- 
lated system, namely layered organic crystals ^-(BEDT- 
TTF)2CsZn(SCN)4 by examining the experimental data 
and the proposed model^S. Disorder arises out of ran- 
dom presence of holes in the charge-ordered layers. The 
low field transport in the insulating phase of these crys- 
tals is characterized by an activated process with m=l. 
To =24 K. The TV characteristics were measured at 
different temperatures and found to follow the power- 
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FIG. 3. Scaling of the bias-dependent conductance data in 
a correlated, charge-ordered organic salt at different temper- 
atures. The original I-V data are taken from Ref. [l^. The 
voltage scale Vo is chosen corresponding to E(Vo) = 2Eo. The 
solid line is a fit to the scaled version of Glazman-Matveev 
expression^: E/E^ = 1 + 0.97q-^^^ + 0.0025?^"/^ + 1.5 x 
lQ-»q™/^ (g = V/Vo). Inset shows two log-log plots of Vo 
vs. Eo (closed symbols) and T (open symbols). The solid 
line is a linear fit to the data (closed symbols) with the slope 
as shown. 



law at large bias with a large exponent (r^ 8.4) simi- 
lar to observations in many systems including conduct- 
ing polymers^, a-Ge^i or a-Si:H2^. Authors suggested a 
model based on unbinding of electron-hole pairs by ther- 
mal excitation in the background of the charge-ordered 
ground state. The power-law characteristics was given by 
/ ^ "l/i +Ua/2kT ^j^g^^ ^ Uo/2kT. Furthermore, the 
onset (or, crossover) field was given by the temperature- 
independent Fo ~ Uo/eX implying xt = 0. Here Uo is the 
strength of the Coulomb potential and A is the screening 
length. Thus, according to the model xt I/^t- Exper- 
imental data adapted froni22, and normalized according 
to Eq. (H} arc presented in Fig. 3. The high quality of 
the collapse of the data at various temperatures points 
toward validity of the scaling to this correlated system. 
Here e was 1 for determining Vo- Ohmic conductances 
at lower T were obtained from extrapolation of the data 
at higher T. The scaled data in the figure represent es- 
sentially the scaling function The inset shows plots 
of Vo vs. Eo and T. The straight line through the closed 
circles proves the validity of Eq. ^ (over a conductance 
range of 14 decades!) and yields a non-zero exponent 
Xt = 0.14 ± 0.01 in contradiction with the model pre- 
diction. Also shown in the figure is a dashed line par- 
allel to the scaled curve at high bias and having a slope 
of 1/0.125 close to 1/0.14 as expected according to Eq. 

A was determined using Fo- The inset in Fig. 3 
shows that Fq increases with temperature. Consequently, 
A must decrease with temperature. Therefore, this casts 




FIG. 4. The nonlinearity exponents xm in several (a) un- 
correlated and (6) correlated classes of materials. The ex- 
ponents are quantized and given by xm ~ 0.08n. This is 
highlighted by plotting xm/0-08 against n with lines dropped 
off the data points to y-axes. Each panel shows a solid line 
passing through the origin with a slope of unity. The inset in 
panel (a) displays the histogram of the exponents in the same 
panel. 



doubt on the conclusion ( i.e., the long-ranged Coulomb 
interaction) drawn from a particular value of Fo- The 
model expression of a temperature-dependent zt is un- 
tenable, particularly at low T or/and at high fields when 
the conductivity becomes temperature-independent i.e., 
'activationless'. If hopping is considered then the con- 
ductivity is given hy a ~ exp(— 2r^/^) where is the 
hopping length and ^ is the localization length. Using 
~ Uo/eF corresponding to the logarithmic potential 
one gets In cr ~ —2Uo/eF£^, an exponential rather than 
power-law dependence. Similar predicament also arises 
in VRH^. We note that for a slightly different poten- 
tial well U = Uo[l — exp(— r/A)] the energy minimization 
leads to r = Xln{Uo/eFX) yielding a power-law depen- 
dence, cr « [F/Foy^ where Fo = Uo/e\ and zt = 2A/C. 

Comparison of the insets in Figs. 1 and 3 highlights 
the difference in behavior of Vo{T) in weak and strong 
localization respectively. Similarly, the huge difference 
in the ranges of y-axes of Figs. 1 and 3 illustrates the 
effects of XT being greater and less than 1. A puzzling 
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fact however is that the scaled curve could be nicely 
fitted by a modified Glazman-Matveev expression^ 
requiring only three powers. 



III. DISCUSSION 

Values of the nonlincarity exponent extracted in this 
work from various 3D strongly localized systems reported 
in literature are given in Table I (Appendix B). Expo- 
nents in conducting polymers and manganites already re- 
ported in Refsj^i^ are shown for convenience in Tables II 
and III. It is quite apparent from the tables that in stark 
contrast to thermodynamic critical phenomena, the uni- 
versal scaling lacks an universal exponent. Recall that in 
2D weak localization oc 1/ i?n . In 3D strong localiza- 
tion one notices the large spread in values of the nonlin- 
carity exponent not only in different systems but also in 
the same material. For example, the five different sam- 
ples of a-Ge examined had for the exponent two values of 
0.16 and three values of ~0.24 (Table I). No correlation 
seems to exist between xt and To(10''K): (0.16, 8), (0.16, 
8.9), (0.24, 0.5), (0.23, 2.6), and (0.24, 1) respectively In 
doped polypyrrole, there are exponents with even differ- 
ent signs, i.e., 0.23 and -0.33. However, all values so far 
lie within the bounds, -1 and -1-1. On a closer inspection 
of the exponents in the Tables we observe a remarkable 
fact that an exponent, positive or negative, could be de- 
scribed as an integer multiple of a number « 0.08 i.e., 
xm ~ 0.08n where n is an integer which can take both 
positive and negative values including zero. The num- 
ber, 0.08 is thus the largest common factor among the 
exponents in the tables. To highlight the quantized na- 
ture of the exponents, xa//0.08 is plotted against n for 
weakly- or non-correlated materials in Fig. 4a, and for 
correlated materials in Fig. 4b with lines dropped off 
the data symbols to y-axes. The solid line in each figure 
is a line with a slope of unity and passing through the 
origin. Points arc seen to be lying close to the lines con- 
firming the quantized behavior and distinct bands form 
around integer values on the y-axis. The width and den- 
sity of a band reflect the spread in the values, and the 
frequency, respectively of the associated exponent found 
in the systems examined so far. The fact that the ex- 
ponents follow the same relation in both correlated and 
uncorrclated systems makes an impressive statement on 
the new universal physics at play in the field-dependent 
conduction in strongly disordered systems. 

Some points are in order, i) The inset in Fig. 4a shows 
a histogram of the exponents in the panel (a). The distri- 
bution exhibits a pronounced peak at n=2 and is clearly 
asymmetric around the peak with a tail towards large 
71. ii) Both conducting polymers and amorphous/doped 
semiconductors follow VRH-type transport. Yet, the ex- 
ponents in the former assume both positive and nega- 
tive values whereas those in the latter assume only pos- 
itive values; iii) For xm = 0, the field-dependence and 




FIG. 5. Summary of the nonlinearity exponent in various 
regimes in the r=0 scaling diagram with same axes- . The 
solid point representing the metal-insulator transition (MIT) 
is predicted to have an another one in xt- The star (cor- 
responding to Ra « 30kn) represents the zone that sepa- 
rates the metallic and insulating states on the 2D curve. The 
dashed line schematically denotes cross-over from 2D to 3D 
with increase of thickness. '?' denotes experimentally either 
unavailable or unverified information. 



temperature-dependence of conductance becomes com- 
pletely separate (Eq. [1]); iv) The field scales remain un- 
affected by an applied magnetic field i.e. xb = (Tables 
I and II). Although a magnetic field does not impart any 
energy to electrons the zero value of xg is still significant 
because a magnetic field destroys time reversal invari- 
ance and is known to influence a wavefunction. However, 
in 2D weak localized regime Eq (O predicts a negative 
exponent, xb = —0.5. There is yet no experimental con- 
firmation of this prediction ; v) For disorder as a control 
parameter, xd in polpyrrole was found to be negative, 
-0.33. Interestingly, Eq (O does yield a negative expo- 
nent, XD = —0.5. On the other hand, Larkin et al.— 
discussed a model of the MIT regime that yields a posi- 
tive exponent, xd = ^ + V where the exponent rj relates 
the energy scale with length scale; vi) Limited data {xd, 
Xb) indicate that xm follows the same quantized behav- 
ior even when M is not temperature. 

Fig. 5 summarizes the experimental findings about the 
nonlinearity exponent in disordered systems in 2D and 
3D on the r=0 scaling diagram^. For highly conduct- 
ing 3D samples, the nonlinear exponent is undefined. As 
the MIT is approached from the metallic side, xt > 1— . 
As MIT is just crossed to the insulating side, Xt < 1— . 
Therefore, the MIT is also the point where xt appears to 
undergo a sharp transition in value from one > 1 to one 
< 1. However confirmation about the transition must 
await further experimental as well as theoretical input. 
For large conductance in 2D, the exponent xt 9. 
Information about xt at larger disorder is lacking. But 
if it is found to have values < 1 as in 3D, an intrigu- 
ing question arises: does the point (star) on the 2D 
curve corresponding to Ru « 30kri represent simply a 
crossover zone or a transition manifest only in the non- 
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linear regime? The dashed hne represents schematically 
a possible path corresponding to increasing thickness of 
a 2D film when the exponent decreases from a large value 
to one closer to 1 as discussed earlier (Fig. 2). 

In summary, heuristic arguments were presented for a 
universal one parameter scaling of field-dependent con- 
ductance in 2D and 3D disordered systems, supported 
by abundant experimental evidence. The associated ex- 
ponent possesses a spectrum of values characteristic of 
various disordered regimes. The unexpected quantiza- 
tion of the exponent in strong localization points toward 
a new physics and calls for further theoretical efforts. A 
possible transition in 2D is speculated and several unre- 



solved issues are discussed. It is apparent that the electric 
field along with the traditional magnetic field should be 
valuable probes for transport in disordered systems. 
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TABLE I. The nonlinearity exponents xm for different materials (mostly exhibiting variable range hopping) obtained using 
the scaling formalism. T,D and B represent temperature, disorder and magnetic field respectively as the control parameter. 
Data sources are shown, m is the VRH exponent and n is the quantization integer.) 



Systems 


m 


XT Xb 


n 


Data Source 


a-Semicond. 












a-Si:Y 


0.5 


0.063 


1 


F Ladicu ct al., Phys Rev B 61, 8108 


(2000) 


a-Si:H 




0.065 


1 


C E Ncbcl ct al., Phys Rev B 46, 6803 


(1992) 


a-SiO 


0.25 


0.072 


1 


A Servini ct al.,Tliin Sol Films 3, 341 


(1969) 


a-CN:H 


0.25 


0.16 


2 


S Kumar et al., J Non-Cryst Solids 338-340,349 


(2004) 


a-Ge 


0.25 


0.155 


2 


M Morgan et al., Phil Mag, 23, 661 


(1971) 


a-Ge 


0.25 


0.16 


2 


M Telnic et al., Phys Stat Sohdi(b) 59, 699 


(1973) 


a-Ge 


0.25 


0.244 


3 


P J Elliot et al., AIP Conf. Proc. 20, 311 


(1974) 


a-Ge 


0.25 


0.23 


3 


T Suzuki ct al., J Non-Cryst Solids 23, 1 


(1977) 


a-Ge 


0.25 


0.24 


3 


N Croitoru et al. Thin Sol Films 3, 269 


(1969) 


a-Gc:Cu 


0.5 


0.37 


5 


A N Aleshin et al., Sov Phys Semicond. 21, 466 


(1987) 


Doped-semicond. 












n-Si:As 


0.5 


0.081 


1 


C Gang et al.. Solid State Comm 72, 173 


(1989) 


n-Si:P 


0.25 


0.15 


2 


T N Rosenbaum et al.. Solid State Comm 35, 663 


(1980) 


n-Si:Mn 


0.5 


0.15 


2 


A V Dvurechcnskii ct al., JETP Lett. 48, 155 


(1988) 


n-GaAs 


0.5 


0.17 


2 


F Tremblay ct al., Phys Rev B 40, 3387 


(1989) 


n-CdSe 


0.5 


0.16 


2 


D Yu et al., Phys Rev Lett 92, 216802 


(2004) 


n-Si:As 


0.25 


0.30 


4 


R W van der Heijden et al., Phil Mag B 65, 849 


(1992) 


p-Si:B 


0.25 


0.22 


3 Y 


M Shwarts et al., Sem Phy, Quan Elec Opto 3, 400 


(2000) 


n-GaAs 


0.5 


0.31 


4 


D Redfield, Adv Phys 24, 463 


(1975) 


n-ZnSc 


0.5 


0.31 


4 


I N Timchenko ct al., Sov Phys Semicond. 23, 240 


(1989) 


d-Gc:Ga 


0.5 


0.5 


6 


T W Kenny et al., Phys Rev B 39, 8476 


(1989) 


n-Si:P 


0.5 


0.64 


8 


J Zhang et al., Phys Rev B 57, 4472 


(1998) 


Composite 












SGNT-PMMA 


0.5 


0.39 


5 


J M Benoit et al., Phys Rev B 65, 241405 


(2002) 


C-PVC 


0.67 


0.63 


8 


L J Adriaanse et al., Synth Metals 84, 871 


(1997) 


Double perovskite 












Ba2MnReO(i 


0.5 


0.155 


2 


B Fisher et al., J Appl Phys 104, 033716 


(2008) 


Organic crystal 












6I-(BEDT-TTF)2 


1 


0.14 


2 


Takahide et al., Phys Rev Lett 96, 136602 


(2006) 


CsZn(SCN)4 












6I-(BEDT-TTF)2 


1 








Takahide et al., Phys Rev Lett 96, 136602 


(2006) 


CsZn(SCN)4 













TABLE II. The nonlinearity exponents for different conducting polymers from Phys. Rev. B 84, 054205(2011). 



System 


m 


XT 


Xd 


Xb 


n 


PPy (powder) 


0.25 


-0.329±0.014 






-4 


PPy (film) 


0.33 


-0.155±0.012 






-2 


PPy (film) 


0.5 


-0.228±0.008 






3 


PPy (film) 






-0.3 




4 


PPy (film) 














PEDOT 


0.5 


-.158±0.003 






2 


PDA 




-0.51±0.02 






6 


PA(nanofiber) 




-0.33±0.02 






4 
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TABLE III. The nonlinearity exponents for OMR phases in different maanganites from Phys. Rev. B 86, 165104 (2012). PI 
and FM represents the paramagnetic insulating and ferromagnetic phases. 



System 










Smo.55(Sro.5Cao.5)o.45Mn03(sc) 












Smo.55Sro.3375Cao.ii25Mn03 


0.17 ±0.01 


2 


0.09 ± 0.01 


1 


Smo.56Sro.45Mn03 


-0.23 ± 0.01 


-3 


-0.14 ±0.01 


-2 


Lao.275Pro.35Cao.375Mn03 


-0.09 ± 0.01 


-1 








Lao.87(Mn203)().i3Mn03 


-0.83 ±0.01 


-10 


-0.16 ±0.01 


-2 


Lao.75Cao.25Mn03 


0.27 ± 0.04 


3 


0.27 ±0.04 


3 


Lao.76Cao.26Mn03/BaTi03 


-0.70 ± 0.01 


-9 


-0.15 ±0.03 


-2 



